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OPTIMUM  LAUNCH  VEHICLE  SIZING 


II 


I.  INTRODUCTION 

In  launch  vehicle  sizing.  It  is  often  necessary  to  determine  the  relation- 
ships between  the  parameters  of  vehicle  size  and  the  parameters  of  mission 
requirements  for  multi-stage  rockets  in  powered  ascent.  Relationships 
of  this  type  are  generally  arrived  at  as  a result  of  an  optimization 
procedure  based  upon  some  reasonable  criterion,  such  as  maximizing  the 
payload  weight  to  launch  weight  ratio,  or  minimizing  the  total  launch  vehicle 
cost.  In  a particular  situation,  the  mission  may  be  to  attain  certain 
launch  altitude  and  velocity.  Then,  if  a set  of  preliminary  design 
parameters  are  also  known  (such  as  the  number  of  stages,  the  specific  impulse 
of  the  stages  and  the  structure  ratios  of  each  stage),  the  output  of  the 
optimization  program  will  provide  sizing  parameters  for  each  of  the  launch 
vehicle  stages. 

The  major  forces  experienced  by  a multi-stage  rocket  in  powered  ascent  are: 
a thrust  force  from  the  rocket  motors,  aerodynamic  forces  of  drag  and 
lift  and  a gravitational  force  of  attraction  towards  the  earth.  Drag  and 
lift  forces  are  nonlinear  functions  of  vehicle  velocity,  while  the  gravitation 
force  depends  upon  vehicle  altitude  and  is  also  a nonlinear  function  of  altitude. 
Because  of  these  facts,  a general  analytical  optimization  procedure  for 
launch  vehicle  sizing  taking  Into  account  the  aerodynamic  forces  and  the 
gravitational  force,  remains  an  unsolved  problem.  For  cases  where  precise 
values  of  vehicle  sizing  are  required,  designers  have  resorted  to  large 
scale  numerical  computation  programs  which  examines  the  value  of  each  para- 
meter Individually  while  fixing  the  value  of  remaining  parameters.  While 
this  numerical  method  does  provide  more  accurate  values,  compared  to  those 
from  an  optimization  procedure  to  be  discussed.  It  Is  time  consuming,  tedious 
and  does  not  provide  a picture  of  the  relationships  between  the  parameters 
of  launch  vehicle  size  and  the  parameters  of  mission  requirement.  In  addition. 

In  a preliminary  vehicle  design  effort  an  optimization  approach  based  upon 
a simplified  model  often  provides  results  of  sufficient  accuracy. 


In  this  report,  a brief  description  of  an  optimization  procedure  for  launch 
vehicle  sizing  for  multi-stage  rockets  In  powered  flight  based  upon  a 

Note:  Manuscript  submitted  February  12, 1979. 
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simplified  model  will  be  presented.  This  will  be  followed  by  an  outline 
of  the  analytical  relationships  for  a Fortran  program  which  can  be  used 
to  compute  preliminary  launch  vehicle  parameters. 

II.  OPTIMIZATION  OF  PAYLOAO  WEIGHT  TO  TOTAL  LAUNCH  HEIGHT  RATIO 


In  this  section,  we  briefly  describe  an  optimization  procedure  due  to 
Ruppe1  which  can  be  used  to  compute  the  launch  vehicle  parameters  for  the 
case  when  the  payload-to-launch  weight  ratio  Is  optimized  (maximized).  The 
analytical  procedure  Is  developed  by  assuming  the  step  rocket  Is  In  vertical 
ascent.  In  vacuum,  under  constant  thrust  per  stage.  Under  these  assump- 
tions, the  gravitational  and  aerodynamic  forces  are  neglected  so  that  the 
assumed  model  Is  not  realistic.  However,  for  Initial  vehicle  sizing,  the 
results  obtained  do  provide  useful  Information. 


The  optimum  weight  distribution  for  a muitl-stage  or  step  rocket  In  vertical 
ascent,  under  constant  thrust  (per  stage).  In  vacuum  will  be  outlined.  First 
the  following  definitions  are  required: 

wQ(1)  * weight  of  vehicle  at  Ignition  of  step  No.  1 
wc(1)  « weight  of  vehicle  at  ctttoff  of  step  No.  1 
Wl(1)  - payload  weight  of  step  No.  1 

* weight  of  vehicle  at  Ignition  of  step  No.  1 + 1 


Wj(1)  ■ weight  of  structure  of  step  No.  1 
Figure  1 shov/s  a 3-stage  multi-stage  rocket. 


In  a multi-stage  rocket,  the  first  stage  accelerates  to  a burnout  velocity 
, at  which  time  the  exhausted  first  stage  structure  Is  ejected  and  the 
second  stage  Ignition  starts.  The  second  stage,  which  starts  with  Initial 
velocity  Vj,  accelerates  the  vehicle  by  a velocity  Increment  of  so 
that  at  the  end  of  stage  two  burn,  the  vehicle  velocity  Is  Vj  + V^.  This 
process  continues  for  additional  stages  until  the  final  burnout  velocity  of 
the  payload  Is  given  by  Vj  ♦ V2  + . . . VN,  where  N Is  total  number  of 
stages.  To  continue,  the  following  definition  will  also  be  required: 

MASS  ratio  of  step  No.  1 
W (1) 

r,  «-S_  (1) 

1 wc(i) 
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3 Stag*  Rocket 


Payload  ratio  of  step  No.  1 

W,(1) 

XJ  = -! 

' «,«> 

Structure  ratio  of  step  No.  1 
Ws(1) 

C - 5 


(2) 


w0(D 


(3) 


From  the  definition 

W1  (N)  = XnH0(N)  = XN  VlWo(N_1) 
N 

= n x4  w ( i ) 


l 


1 O' 


(4) 


N = Total  number  of  stages 

The  ratio  of  total  launch  weight  to  payload  weight  can  now  be  written  as 
1 


(5) 


V1) „G 

Wl(N)  xlx2  * * * * XN 

Where  G Is  called  the  growth  factor. 


The  objective  of  the  optimization  procedure  Is  to  minimize  the  growth 
factor,  G,  so  as  to  minimize  the  launch  vehicle  weight  to  payload  weight 
ratio.  This  Is  equivalent  to  maximizing  the  payload  weight  to  launch 
weight  ratio. 

In  the  Ruppe  optimization  procedure,  certain  mission  parameters  are 
assumed  known.  These  Include  the  Ideal  velocity  (Vjd),  specific  Impulse 
of  each  stage  (SPI),  structure  ratio  of  each  stage  (E),  and  the  Initial 
load  factors  (A). 
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The  Ideal  velocity  Is  defined  by: 


V, . = V + V + V + V . + V„  + V 
id  m w er  si  9 a 


where  Vm  Is  the  desired  mission  velocity,  Vw  Is  the  velocity  component 

allocated  for  launch  windows  capability,  Ver  Is  the  velocity  component 

due  to  earth  rotation,  Vsl  Is  the  velocity  component  for  system  losses 

(such  as  variations  In  motor  efficiency),  Vg  Is  velocity  component  due  to 

gravity  loss  and  V Is  velocity  component  due  to  aerodynamic  effects,  such 
•» 

as  drag  and  lift. 

Then,  the  optimization  procedure  first  solves  for  the  unknown,  XN  from 


nil  - 


Cj  E..  ^ An 


where  R Is  computed  from 
R * exp  (V1d/c) 


With  C being  the  average  exhaust  velocity  of  the  stages  given  by 
_ i 


where  Cj  Is  the  exhaust  velocity  of  gas  of  jth.  stage.  Under  our 
assumption  of  vacuum  and  constant  gravitational  acceleration. 


ci  ’ SPI(i)V 


where  SPI^  Is  the  specific  Impulse  of  jth.  stage  and  gQ  Is  the 
acceleration  of  gravity  at  earth's  surface. 


Since  in  Eq.  7,  X^  Is  the  only  unknown,  the  value  of  XN  must  be  solved 
by  an  Iterative  procedure.  Once  this  is  accomplished,  X^,  1 » 1,  2 . . 
can  be  obtained  from 
E 


1 


Xi  = 


Si  / xn  \ 

Ci  \EN  + XN  ' 


Ei* 


01) 


1 - 


. N-l, 


With  X.  now  known,  the  calculation  can  proceed  to  determine  the  parameters 
of  vehicle  sizing  in  the  following: 

Stage  Weight  (lbs) 


W0(l)  = G-PAY(N) 


(12) 


Fuel  Weight  (lbs) 

Wf ( 1 ) = wo(i) 


f1  = 


xi  + Ei 


Structure  Weight  (lbs) 


Ws(1)  - W0(1)  * E1 


Payload  Weight  of  1th  stage  (lbs) 
PAY ( 1 ) *=  W0(1)  - Wf(1)  - ws(i) 

Exhaust  Velocities  (kilome ter/sec) 


SPI 


(D 


(13) 


(14) 


(15) 


(16) 
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Thrust  (lbs) 


TO)  = A(1)Wo(1)  (17) 

Burn  time  (sec) 

, t Wf(1) 

TbO)  - -£ * SPIm  (18) 

b T(1) 

III.  FORTRAN  PROGRAM 

Rally  I (Section  V)  Is  a Fortran  program  which  computes  the  launch 
parameters  of  a step  rocket  for  the  case  when  the  payload  weight  to 
launch  vehicle  weight  Is  maximized.  The  program  Is  based  upon  the  model 
which  assumes  the  step  rocket  Is  In  vertical  ascent.  In  vacuum,  and  under 
constant  thrust  per  stage.  Although  these  assumptions  may  seem  somewhat 
restrictive,  the  results  obtained  do  provide  first  order  approximations 
for  launch  vehicle  design. 

A list  of  the  program  INPUTS  and  OUTPUTS  will  now  be  provided: 


INPUTS 

A.  Structure  ratios,  E^,  entered  as  data  In  step  200  under  EPS 
for  values  corresponding  to  stages  1,  2,  3,  ...  N,  In  order. 

B.  Staqp  exhaust  velocities  (kilometer/sec),  entered  as  DATA  In 
step  200  under  C for  stages  1,  2,  3,  . . . N In  order. 

C.  Specific  Impulses,  entered  as  DATA  In  step  210  under  SPI_  for  stages 
1,  2,  3,  . . . N In  order. 

D.  Initial  load  factors  (ratio  of  thrust  to  stage  weight)  entered  as 
DATA  In  step  220  under  AK  for  stages  1 , 2,  3,  ...  N In  order. 

E.  Payload  weight,  entered  In  step  400  under  PAY  In  lbs. 

F.  Number  of  stages,  entered  In  step  500  under  N.  There  Is  no  restric- 
tion on  the  value  of  N,  except  the  values  In  the  DIMENSION  size 
(steps  100,  110,  120)  must  be  made  equal  to  or  greater  than  N. 
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6.  Ideal  velocity,  (kilometers/sec),  entered  In  step  600  under  VID. 
OUTPUTS 

The  OUTPUTS  of  RALLY  1 are: 

A.  The  growth  factor,  6 

B.  Thrust  (lbs)  per  stage 

C.  Stage  weight  (lbs)  per  stage 

D.  Fuel  weight  (lbs)  per  stage 

E.  Payload  weight  (lbs)  per  stage 

F.  Structure  weight  (lbs)  per  stage 

6.  Specific  Impulse  (seconds)  per  stage 
H.  Burn  time  (seconds)  per  stage 
IV.  NUMERICAL  EXAMPLES 

In  this  section,  a sample  calculation  will  be  given  to  show  the  use  of 
FORTRAN  Program  RALLY  1.  For  convenience,  we  consider  a 2-stage  rocket 
whose  mission  Is  to  launch  a payload  Into  a circular  orbit  at  an  altitude 
of  100  nautical  miles.  The  mission  parameters  for  this  example  are 
given  In  a U.  S.  Air  Force  document  titled:  Space  P1anrter*s  Gulde^,  (SPG), 
page  V-20.  In  this  document,  a graphic  procedure  for  vehicle  sizing, 
which  maximizes  the  payload  weight  to  launch  vehicle  weight.  Is  given. 

The  results  of  this  graphic  procedure  will  be  compared  to  the  values 
obtained  from  Rally  1. 

A.  Mission  Parameter 


The  following  mission  parameters  were  given  In  SPG 


Orbit 

100  Nautical  Miles,  Circular 

Number  of  stages 

2 

Payload  Weight 

81,324  lbs 

Ideal  velocity 

32.480  ft/sec 

Vacuum  specific  impulse 

Staqe  1 
395  sec 

Staqe  2 

419  sec  ; 

Thrust  to  weight  ratio 

173 

1.0 

TABLE  1 - Mission  Parameters  (SPG) 


E 
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B.  Rally  1 Parameters 


In  order  to  use  Rally  1,  the  mission  parameters  of  Table  1 are  converted 
Into  appropriate  units  shown  In  Table  2. 


Number  of  stages 


Payload  Weight 


Ideal  Velocity 


Exhaust  velocity 
Load  factors 
Structure  ratio 
Specific  impulse 


81,324  lbs 


9 .885  km/sec 


STAGE  1 


3.871  km/sec 
1.3 
.0665 
395 


STAGE  2 


4,1062  km/sec 
1.0 
.0965 
419 


TABLE  2 - Rally  1 Inputs 


In  Table  2,  the  Ideal  velocity  was  obtained  simply  by  conversion  of  units 
Into  km/sec.  The  exhaust  velocity  was  obtained  by  using  Eq.  16. 

Cj  . SPI(j)  • g0  (18) 

where  Cj  Is  the  exhaust  velocity  of  the  jth  stage,  gQ  Is  the  acceleration 
of  gravity  at  earth's  surface  (km/sec),  and  SPI^j  Is  the  vacuum  specific 
Impulse.  The  values  of  "structure  ratios",  E,  are  computed  from  the  values 
of  "structure  factor"  and  the  values  of  stage  weights  shown  In  Table  3 for 
the  SPG  case.  Since  the  definition  of  "structure  factor"  Is 
W (1) 

Sf(1)  - -5 (19) 

where  S^(1)  Is  the  structure  factor  of  1th  stage,  W ( 1 ) Is  the  structure 
weight  of  1th  stage  and  W'(1)  Is  the  stage  weight  of  1th  stage  only  (does 
not  Include  weight  of  stages  1 + 1 , 1 + 2,  . . . N);  the  relationship 
between  E(1)  and  S^(1)  can  be  determined  by  noting 


PARAMETERS 


Stage  v/eight  (lbs) 
Structure  factor 
Structure  Ratio 
Propellant  Weight 
Load  Factor 
Thrust  (lbs) 

Stage  Velocity  Inc: 


Growth  factor 


SPG 

STAGE  1 

Staqe  2 

1 .854m 

,386,000 

0.0840 

.0965 

1.340  m 

275JOOO 

1.3 

1.0 

2.41m 

284  000 

16,  050 

16,400 

22.798 


RALLY  1 


Stage  llStaqe  2 


1.875m  1360,045 


.0965 

1,389m  | 244,343 

> 1-° 
2.44m  360,450 

17,197  15,285 


23.057 


Table  3 - Comparison  of  SPG  and  Rally  Results 


W0(1)  » W'(l)  + W-(2) 

Since 

Wq(2)  - W'(2) , W'(1)  = W0(l)  - W-(2)  (20) 

The  structure  factor  of  Eq.  19  becomes 


sf(l) 


ws(l) 

Xo0)  - W'(2) 


(21) 


Dividing  by  WQ(1),  we  note  the  relationship  between  E and 

/ Mn(2>  \ 

E(l)  - Sf(l)  (l  - -2 ) (22) 

f V w0(i)  / 

E(2)  - Sf(2)  (23) 

Thus  E (1)  and  E (2)  can  be  computed  from  the  values  of  SPG  data  given 
In  Table  3. 

Sf(l)  - .0840 

Sf(2)  - .0965 

WQ(1)  - 1,854,000  lbs 

WQ(2)  = 386,000  lbs 

/ 

The  computed  values  of  launch  vehicle  parameters  using  Fortran  Program 
Rally  1 are  shown  In  Table  4.  A summary  of  the  computations  are  shown  In 
Table  3 for  comparison,  and  reasonable  agreement  between  the  SPG  values  and 
Rally  1 values  ca(.  be  seen  for  the  2-stage  rocket  under  consideration. 
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OBAR RBAR 

• • -3.98860 • • -1-1. -921-14 


X LAMM 

0.-70042  0.22562 


LAMPDA  MUMPER  *1-  -IS 0.-19223 

LAMBDA  MUMPER  2 IS 0. -22562 


L-DAD  -RAT- ID  -MUMPER  -1  -IS 3 .-865 01 

LORD  -RATIO  -MUMPER  -2  IS 3.10446 


DELTA  -V  MUMPER  I -IS  • • • -5.-23345 
DELTA  -V  -MUMPER  2 -IS 4 .-65 1-66 


GROWTH  -FACTOR  -IS  • • -23 . 05693 


PARAMETERS  FOR  STAGE  1 FOLLOW 

THE  -THRUST  -IS  • -2437606.7  PDUNDS 
THE  -STAGE  WEIGHT  IS  • -1875082.1  -POUNDS 
THE  -FUEL-  WEIGHT  -IS  • -1389938. -6  -PDUMDS 
THE  -PAVLDAD  -WEIGHT-  -IS  • • -360450.5  -POUNDS 
THE  -STRUCTURE  -WEIGHT  -IS  • • -124693.  0 -POUNDS 
THE  -SPECIFIC  -IMPULSE  -IS  • -395.  00000  -SECONDS 
THE  -BURN  -TIME  -IS  • -225.231-47  -SECONDS 


PARAMETERS  -FOR  STAGE  -2  -FOLLOW 


THE  -THRUST  -IS  • • -360450.5  -POUNDS 

THE  -STAGE  WEIGHT  -IS  - - -360450.5  POUNDS 

THE  -FUEL-  -WEIGHT  4S  • • -244343.  1 POUNDS 

THE  -PAYLOAD  -WEIGHT  -IS 81324.  0 -PDUNDS 

THE  -STRUCTURE  -WEIGHT-  -IS 34783.  5 -PDUMDS 

THE  -SPECIFIC  -IMPULSE  -IS  • -419.00000  -SECONDS 
THE  -BURN  -TIME  -IS  • -284.  03271-  -SECONDS 


TABLE  4 - Rally  1 OUTPUTS 
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sssssssssr^ 
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I 


04600 



•••DO  50  1-lfM 

04700 



• • -ALAM  <I>  ■ CEPS  CD>  <1 .--  <0  CM) vC  ci»  * CAL-AMN/1  C CEPS  CM>  .fiLftMM>  > > > > 

) . 



* • -1* 

04800 



• • TYPE  5rIrALAM<I* 

04900 

..•■I** 

• • -FORMAT  C'-  LAMBDA  -NUMBER  'rWr'-  -IS  -'-rF10.5> 

05000 

• • *50  • 

• • CONTINUE 

051  00 



• • TYPE  9 

05200 



• • DO  -60  I=l-rN 

05300 



• • RATIOCI)— l./*CALAM<I>-t-EPS<I>* 

05400 



‘ - TYPE  -6  r I r RAT  1 0 < I*  i 

05500 

• • FORMAT  C'-  LOAD  RATIO  MUMBER  -'-rllr'-  IS  -'-rF10.-5>- 

05600 

• -60  • 

• • -CONTINUE 

05700 



• • TYPE  -9 

05800 



• • DO  -70  -I=I-»-N 

05900 



• • -DELV  C I*=0  C D-.ALD6  CRAT  ID  CI» 

06000 



• • TYPE  -7r  IfDEL-V  CI> 

06100 

• • •?  • • 

• • -FORMAT  C'-  DELTA  V NUMBER  'rU  r'-  -IS  -'-rF10.-5>- 

06200 

• • ?0  • 

• • -CONT I NUE 

06300 



• • TYPE  -9 

06400 



• • -ALAMT-1  .-0 

06500 



• • DO  -80  1=1  » M 

06600 



• • •ALAMT=ALAMT^ALp.M  CI> 

06700 

. . 80  • 

• • -CDNT I NUE 

1168OO 



• • •G=I->ALAMT 

06900 

...... 

• • -AL-W=G*PAY 

07000 



• • TYPE  8>G 

071  00 

• • *8  • • 

• • -FORMAT  C'-  -GRDWTH  -FACTOR  -IS  -'->F10.-5>- 

07200 

. . .9  . . 

• • -FORMAT  <yy> 

07300 



• • -SW<1*»ALW 

07400 



• • -DO  -90  I-l  rN 

07500 



• • >WF  CI*»SW  <!*♦  < CRAT  ID  < I>~  1*>RAT  I-D  C I» 

076-00 



-••WS<I*=SW<I*.-EPS<I>- 

07700 



•••SW<I--*-l*=SW<I*-CWFCI>--*-WSCI>* 

07800 



• • -THRCD-AKCD-^SWCl* 

079.00 



• • •BTCI*=WF<I>/'THR«I**SPI-<I> 

08000 

. . .90  • 

• • -CONTINUE 

08100 



■ • -TYPE  -9 

08200 



• • -DO  100  I-lfN 

08300 



• • -TYPE  -l-40r I 

08400 

• • 140 

• • -FORMAT-  C'-  -PARAMETERS  -FOR  -STAGE  -'-rllr'-  -FOLLOW '•'*>• 

08500 



• -TYPE  -lSOrT-HRCI* 

08600 

■C 

• '-TYPE  160.SWCI) 

08700 



i . 

• • -TYPE  -1  70rWF  Cl* 

09800 



• • -TYPE  -180rSW  <1+1*  • * 

08900 



• • -TYPE  190  rWSCI> 

09000 



• • -TYPE  200rSPI  CI> 

09100 



•• -TYPE  -21- Or  BTC  I* 

09110 



• • -TYPE  9 

09200 

• • 150 

• • -FDRMAT-  <'•  -THE  -THRUST  -IS  -'-rFl  O.-l-r -POUNDS'* 

09300 

• • *1*0 

• • -FDRMAT  C'-  -THE  -STAGE  -WEIGHT  -IS  -'-rFl  0.1  r'-  -POUNDS') 

09400 

• • 170 

• • -FORMAT-  <'•  -THE  -FUEL  -WEIGHT  -IS  -'-rFlO.-l-r'-  -POUNDS'* 

09500 

• • *1-80 

• • -FORMAT  C'-  -THE  -PAYLOAD  -WEIGHT  -IS  -'-rFl  0.1  r '•  -POUNDS'* 

09600 

. . -190 

• • -FORMAT  C'-  -THE  -STRUCTURE  -WEIGHT  -IS  -'-rFl  0.1  »•'•  -POUNDS'* 

09700 

• • 200 

• • -FORMAT <'•  -THE  -SPECIFIC  -IMPULSE  -IS  -'-rFl-O.-Sr'-  -SECONDS'* 

09000 

• • -210 

• • -FORMAT  </•  -THE -BURN  -TIME  -IS  -'rFl-O.-Sr'-  -SECONDS'* 

09900 

• • -*00 

• • -CONTINUE 

10000 

• • -END 
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